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Abstract. We study the asymptotic behaviour of a class of self- attracting motions on M.'^. 
We prove the decrease of the free energy related to the system and mix it together with 
CO ' stochastic approximation methods. We finally obtain the (limit-quotient) ergodicity of the 

self-attracting diffusion with a speed of convergence. 
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1. Introduction 

1.1. Statement of the problem. This text is devoted to study the asymptotic behaviour 
of a Brownian motion, interacting with its own passed trajectory, so-called "self-interacting 
motion". Namely, we fix an interaction potential function ly : M'^ — )■ M, and consider the 
stochastic differential equation 



dXt = V2dBt -{^ I VW{Xt - X,) d^ dt. 



where (-Bt,t > 0) is a standard Brownian motion, with an initial condition of given Xq (with 
the condition of continuity at t = 0). This equation can be rewritten using the normalized 
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occupation measure fit'- 



f^t = - Sx, ds, 



where 6^ is the Dirac measure concentrated at the point x. Using this convention, the equa- 
tion ([T]) becomes 

(2) dXt = V2dBt-VW* fitiXt) dt, 

where * stands for the convolution. 

Note that the equations ([1]), ([2]) clearly have singularities at t = 0, which is the reason why 
sometimes they are considered only after some positive time r > 0. We discuss the existence 
and uniqueness questions for the solution in the appendix. 

Similar problems have already been studied since the 90's, for instance by Durrett and 
Rogers |8], or Benaim, Ledoux and Raimond [2], initially to modelize the evolution of polymers 
or ants. The first time-continuous self-interacting processes have been introduced by Durrett 
and Rogers [8] under the name of "Brownian polymers" . They are solutions to SDEs of the 
form 

(3) dXt = dBt+n f{Xt- X,) ds"] dt 

where (-Bj,t > 0) is a standard Brownian motion and / a given function. We remark that, 
in the latter equation, the drift term is given by the non-normalized measure t^t and not by 
/it as the process we will study here. As the process {Xt,t > 0) evolves in an environment 
changing with its past trajectory, this SDE defines a self-interacting diffusion, which can be 
either self-repelling or self-attracting, depending on the function /. In any dimension, Durrett 
& Rogers obtained that the upper limit of |Xi|/t does not exceed a deterministic constant 
whenever / has a compact support. Nevertheless, very few results are known as soon as the 
interaction is not self-attracting. 

Self-interacting diffusions, with dependence on the (convoled) empirical measure {fit,t > 0), 
have been considered since the work of Benaim, Ledoux & Raimond [2]. A great difference 
between these diffusions and Brownian polymers is that the drift term is divided by t. This 
implies that the long-time away interaction is less important than the near-time interaction 
(the interaction is not "uniform in time" anymore). Benaim et al. have shown in [21 E] that 
the asymptotic behaviour of fit can be related to the analysis of some deterministic dynamical 
flow defined on the space of the Borel probability measures. Afterwards, one can go further 
in this study and give sufficient conditions for the a.s. convergence of the empirical measure. 
It happens that, with a symmetric interaction, fit converges a.s. to a local minimum of 
a nonlinear free energy functional (each local minimum having a positive probability to be 
chosen), this free energy being a Lyapunov function for the deterministic flow. These results 
are valid for a compact manifold. Part of them have recently been generalized to M*^ (see 
[0]) assuming a confinement potential satisfying some conditions — these hypotheses on the 
confinement potential are required since in general the process can be transient, and is thus 
very difficult to analyze. In these works, no rate of convergence is obtained. Most of these 
results are summarized in a recent survey of Pemantle [12j, which also includes self- interacting 
random walks. 

Coming back to the process introduced by Durrett & Rogers, all the results obtained have 
in common that the drift may overcome the noise, so that the randomness of the process is 
"controlled". To illustrate that, let us mention, for the same model of Durrett & Rogers, the 
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case of a repulsive and compactly supported function /, that was conjectured in [8] and has 
been partially solved very recently by Tarres, Toth and Valko |15] : 

Conjecture (Durrett & Rogers [8]). Suppose that / : M — )■ M is an odd function of compact 
support, such that xf{x) > 0. Then, for the process X defined by ([3]), the quotient Xt/t 
converges a.s. to 0. 

In ([T]), the drift term is divided by t, and so it is bounded for a compactly supported 
interaction W. As for the process of the conjecture, the interaction potential is in general 
not strong enough for the process ([1]) to be recurrent, and the behaviour is then very difficult 
to analyze. In particular, it is hard to predict the relative importance of the drift term (in 
competition with the Brownian motion) in the evolution. 

On the other hand, in our case of uniformly convex W, the interaction potential is attractive 
enough for the diffusion (a bit modified) to be comparable to an Ornstein-Uhlenbeck process, 
what gives an access to its ergodic behaviour. 

Another problem, related to the one considered in this paper, is the diffusion corresponding 
to MacKean and Vlasov's PDE. Namely, consider the Markov process defined by the SDE 

(4) dYt = V2dBt-VW*ut(Yt)dt, 

where z/f stands for the law of Yt, and ly is a smooth strictly uniformly convex function. 

The questions of the asymptotic law for Y have been intensively studied these last years, 
by Carrillo, MacCann & Villani [5], Bolley, Guillin & Villani |1], or Cattiaux, Guillin & 
Malrieu [7] for instance. It turns out that, under some assumptions, the laws Ut converge to 
the limit measure u*. This measure is characterized as a fixed point of a map 11 : z/ i— )■ n(z/) 
associating to a measure z/ the probability measure 

n(z/)(dx) := ^e-'^*"("Ma;, 

which is the stationary measure of the process, with vt in the right-hand side of (jl]) replaced 
by V and Z = Z^vs, the normalization constant. 

In particular, Carrillo, MacCann & Villani [5] have shown, using some mass transport tools, 
that the relative free energy corresponding to vt with respect to v* decreases exponentially fast 
to 0. Then Talagrand's inequality allows to compare the relative free energy to the Wasserstein 
distance in case of uniform convexity of the interaction potential W ^ and so they have obtained 
the decrease to of the quadratic Wasserstein distance between z/j and z/*. 

We remark that a huge difference between the preceding Markov process and the (non- 
Markov) self-interacting diffusion is that the asymptotic a-algebra is in general not trivial for 
the non-Markov process. Nevertheless, we will use a similar mass transport method to show 
the convergence of the empirical measure /ij. 

1.2. Main results. Our results are analogous to those of Carrillo et al. |5]: under some 
assumptions imposed on the interaction potential W ^ we show that the empirical measure [it 
almost surely converges to an equilibrium state, which is unique up to translation: 

Theorem 1 (Main result). Suppose, that W G C'^iW'-), and: 

1) spherical symmetric: W{x) = W{\x\); 

2) uniformly convex: denoting by E>'^~^ the {d — 1) — dimensional sphere, 

(5) 3Cw>0: \/xeR'^,yveS'^~\ V^ >C, 



dv^ 



w, 
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3) W has at most a polynomial growth: for some polynomial P, we have 

(6) \/xeW^ \W{x)\ + \VW{x)\ + \\V'^W{x)\\<P{\x\)- 

Then, there exists a unique symmetric density poo • '^'^ -^ '^+! such that almost surely, there 



exists Coo such that 

i- I *— weakly 



A^i := 7 / OxA^ ^ Poo (a; - c^) dx. 

Moreover, there exists a > such that the speed of convergence of fit toward poo(- + Coo) for 
the Wasserstein distance is at least exp{— a '°+\/logt}, where k is the degree of P. 

Remark 1. The assumption [T]) corresponds to the physical assumption of the interaction force 
between two particles being directed along the line joining them, and to the third Newton's law 
(that is the equality between the action and the reaction forces). The symmetry assumption 
cannot be omitted, as shows an example in the appendix. 

Remark 2. We will suppose in the following, without any loss of generality, that P > 1 is of 
degree k > 2 and such that for all x,y E M"^, we have P{\x — y\) < P{\x\)P{\y\). Indeed, we 
choose P(|a;|) = A{1 + |a;|'^), where A is a constant large enough. This will be used in §2.21 

The origin of the following remark will be clear after the discussion in §2.41 

Remark 3. The density poo is the same limit density as in the result of [S], uniquely defined 
(among the centered densities) by the following property: poo is a positive function, propor- 
tional to e-^*^°°. 

We can also consider the same drifted motion in presence of an external potential V. For 
this, the following result is a generalization of Theorem [T] (where we replace Cw by C in the 
notation) : 

Theorem 2. Let X be the solution to the equation 

(7) dXt = V2dBt - (w{Xt) + - f VW{Xt - X,)ds j dt. 

Suppose, that V G C^R"^) and W G ^^(R'^), and: 

1) spherical symmetric: W{x) = W{\x\); 

2) V and W are convex, lim|2;|^oo V{x) = +oo, and either V or W is uniformly convex: 



#1/ 

3C>0: yxEW^,yv eS"^-^ 



dv"^ 



>C or Vx,Vt;, 



dv"^ 



>C; 



3) V and W have at most a polynomial growth: for some polynomial P we have Vx E Mf^ 
(8) \V{x)\ + \W{x)\ + \VV{x)\ + \VW{x)\ + ||VV(a;)|| + \\V^W{x)\\ < P{\x\). 

Then, there exists a unique density poo : IR'^ — ?• IR+, such that almost surely 

i / _ - *— weakly , , . 

f^t = - ox.ds > Poo{x) dx. 



t ./q t-s>oo 



As the proof of the latter Theorem coincides with the proof of Theorem [T] almost identically, 
we do not present it here. It suffices to add V in the arguments below. Moreover, if V is 
symmetric with respect to some point g, then the corresponding density poo is also symmetric 
with respect to the same point q. 
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The proof of Theorem [1] is spht into two parts. Consider a natural "reference point" for a 
measure /i: 

Definition 1. Consider a measure /i on W^, decreasing fast enough for W * n to he defined. 
The center of /i is the point c^ = c{fi) such that VW * f^{Cfj) = 0, or equivalently, the point 
where the convolution W * fi (the potential generated by fi) takes its minimal value. Also, we 
define the centered measure /i'^ as the translation of the measure /i, bringing c^ to the origin: 

(9) ^,'^{A)=^,{A + c,). 

Remark 4. This notion of center had been previously introduced by Raimond in [13]. Indeed, 
to study the linear attracting c?- dimensional case of Brownian polymers, Raimond has defined 
the center and proved that the process remains close to q = c{nt) (and that q converges a.s.). 
A sufficient condition for the existence of the center is that W is convex, and it is unique if 
W is stricty convex. 

The first part of the proof of Theorem [T] consists in proving the convergence of centered 
occupation measures: 

Theorem 3. Under the assumptions of Theorem [H, for some symmetric density function 
Poo '■ IR'^ — )• IR+, we have almost surely 

f, *— weakly , , . 

/it > Poo{x) dx 

t—^oo 

The second is the convergence of centers: 

Theorem 4. Under the assumptions of Theorem [H almost surely the centers Ct := c{pt) 
converge to some (random) limit Cqo ■ 

It is clear that the two latter theorems imply the main result. Let us sketch their proofs. 
1.3. Outline of the proof and physical interpretation. 

1.3.1. Existence and uniqueness. First, a standard remark is Markovianization: the behaviour 
of the pair {Xt, pt) is Markovian. The reader will find it, together with some other standard 
remarks, in §2.1.1[ Unfortunately, the Markov process {Xt, pt) is infinite-dimensional and, in 
general (except for the case of a polynomial interaction W), we do not manage to reduce to 
a finite-dimensional process. So, we do not use this information directly in order to obtain 
interesting properties on pt, because the state space is then too large. 

After this remark, we discuss the global existence and uniqueness for the solutions of ([2]) in 
g2Tl 

1.3.2. Discretization. A next step is discretization: we take a (well-chosen and deterministic) 
sequence of times T„ — )■ oo, with T„ 3> T^+i — Tn ^ 1, and consider the behaviour of the 
measures pt„- As T„ ^ T„_,.i — T„, it is natural to expect (and we will give the corresponding 
statement) that the empirical measures pt on the interval [T„, T„+i] almost do not change and 
thus stay close to fiT„- So, on this interval we can approximate the solution Xt of ([2]) by the 
solution of the same equation with pt = PTu '■ 

dYt = V2dBt -VW * pTAYt)dt, tG [T„,T„+i], 

in other words, by a Brownian motion in a potential W * fiT„ that does not depend on time. 
On the other hand, the series of general term T^+i — T„ increases. So, using Birkhoff Ergodic 
Theorenu, we see that the (normalized) distribution /i[T„,T„+i] of values oi Xt on these intervals 



see for instance [18], chap. XIII 
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becomes (as n increases) close to the equilibrium measures n(/iT„) for a Brownian motion in 
the potential W * fiT„ , where (see §3.ip 

1 



n(/x)(dx) 



Z(/.) 



g-M/*M^)^^^ Z(/i) := / e-^^'^^^Mx. 



But 



so we have 



T 



/^r„4 



T 



/^T„+i 



2n , -iri+l 



— T 



n+1 



and 






/^T„ + 



n(/x, 



Tn+l — Tn 



T 

J- 7 



n+1 



/^[T„,T„+i] 



/^T„ + 



^ra+1 — Tn 



T 



(n(/iTj -/^T„), 



n+1 



-'n+l ~ -'n J-n+l 



(n(/irj -/irj- 



This could motivate us to approximate the behaviour of the measures Ht by trajectories of 
the flow (on the infinite-dimensional space of measures) 

(10) A = i(n(/x)-^), 

or after a logarithmic change of variable 9 = log t, 

(11) /i' = n(/i) — /i. 

In fact, it is not a priori clear that the flow defined by ( TTTl) exists, as the space of measures 
is infinite-dimensional. Though the flow can be shown to be well defined on a subspace of 
exponentially decreasing measures, we prefer to avoid all these problems by working directly 
with the discretization model in §3.11 Nevertheless, this flow serves very well in motivating 
the considered functions and lemmas describing their behaviour, as the discretized procedure 
we have is in fact the Euler method of finding solutions to (fTTl) . 

1.3.3. Physical interpretation: gas re- distribution. Before proceeding further, let us give a 
physical interpretation to the flow (fTTl) . predicting its asymptotic behaviour. Namely, note 
that a Brownian motion drifted by some potential V, 

dXt = V2dBt - VV{Xt)dt, 

can be thought as movement of gas particles under this potential, and the stationary prob- 
ability measure, m = -^e~^dx, is the density with which the gas becomes distributed after 
some time passes. So, in dimension one, a discrete approximation to the flow ( ITTj) can be seen 
as follows. We take a tube, filled with VT-interacting gas, separated in a plenty of very small 
cells (see Fig. [1]). 



.^_ 



.^_ 



.^_ 



./_ 



-/_ 



Figure 1. Gas: phase "separation'' 
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Each unit of time, small parts (of proportion e) of gas in these cells are separated, allowed 
to travel along the tube, and are proposed to equilibrate in the potential generated. This 
part of all the gas being small, its auto-interaction is negligible, thus their new distribution is 
governed by the field V := W * ^, generated by the major part of the particles staying fixed 
to their cells. The small part is then equilibrated to its weight e times n(/i). 



Figure 2. Gas: phase "re-distribution" 

Then, it is separated again by the cells, thus the distribution after such step becomes 

(1 — e)^ + el[{fi) = /i + e{I[{fi) — /i). 

On the other hand, this procedure does not require any work (in the physical sense) to be 
done: the only actions are opening and closing the doors. So, due to the general principle, 
one can expect that the system will tend to its equilibrium. And a tool allowing to show that 
it is the case is the free energy, that we recall in the next paragraph. 

We conclude by noticing that the same physical interpretation can be considered for the 
problem in any dimension d, by placing in W^~^^ two close parallel walls (corresponding to the 
tube in dimension one), and placing the cells along them. 



1.3.4. Free energy functional. A tool allowing to show the convergence of trajectories of (ITTl) 
is the free energy that, due to a general physical principle, should not increase along the 
trajectories as long as we do not do any work. 

Namely, consider an absolutely continuous probability measure /i = fi{x)dx (by an abuse 
of notation, we denote the measure and its density by the same letter). Imagine fi{x) as the 
density of a gas, particles of which implement the Brownian motion v2di?j, as well as interact 
with the potential W{x — y). Then, one defines the free energy of /i as the sum of its "entropy" 
T-L and "potential energy" : 



(12) 



J^ifi) ■.= n{fx) + - / n{x)W{x-y)fi{y)dxdy, 



where the entropy of the measure /i is 



H{fJ,) := / /i(x) log/i(x)dx. 



(13) 



Then, as we have already said, a general physical principle says that, as we are doing no work 
on the system, the free energy should decrease, and the system should tend to its minimum. 
Indeed, the free energy J-" is a Lyapunov function for the flow fITT]) (when it is defined, though 
it is defined only for measures that are absolutely continuous with respect to the Lebesgue 
measure, and otherwise, J-'(/i) = +oo). This can be seen by joining two statements: on one 
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hand, the measure 11 (/i) is (what corresponds to the same physical principle) the unique global 
minimum of a free energy 

jr^(/i) := n{i2) + / V{x)i2{x) dx, 

of a non-interacting Brownian motion in the exterior potential V = W * fi (see §1.3.31 and 
Lemma [7] in §2.4p . The second is the inequality 

where m = n(/i). On one hand, it can be easily seen by an explicit computation, noticing that 
the entropy part is convex. On the other hand, such a differentiation corresponds to replacing 
some small parts of the gas distributed with respect to the measure fi by the one distributed 
with respect to the measure m, and in the right-hand side we have the corresponding free 
energies of these small parts in the potential, generated by the main part of the gas. 

Then, differentiating the function J-" along the trajectories of the flow flTTl) . one finds for the 
solution fi{9) 

— J'(/i(6')) < J'wMe)O^Wn^{e)iM)) - J^wMS)iM) < 0, 

with the equality if and only if fi{9) = n(/i(6')). 

Finally (and we recall these arguments in §3.ip . the fixed points of 11 are exactly the trans- 
lation images of the density poo, that is the centered global minimum of the functional J-'. 
So, roughly speaking, the function T is the Lyapunov function of the flow ( ITTj) . The words 
"roughly speaking" here refer to that these arguments are non-rigorous: we avoided showing 
that the flow is indeed well-defined, the free energy functional is defined only for absolutely 
continuous measures, etc. Though all of this serves well as a motivation to (rigorous) lemmas 
of free energy behaviour used in this paper. 

We conclude this paragraph by indicating that for the dynamics in presence of an exterior 
potential V (the case of Theorem |2]) one has to replace the free energy function by 

J^v,w{fJ') ■= 'H{fi) + V{x)fi{x)dx + - / n{x)W{x-y)n{y)dxdy. 

and, instead of J^w*^J.y consider J^v+w*fi for the energy of "small parts". 

1.3.5. Conclusion. We are now ready to conclude the sketches of the proofs of Theorems [3] 
andm (as it was already mentioned, they immediately imply Theorem [1]). 

Namely, we consider the discretized Euler-like evolution of the flow ( TTOl) . defined by the rule 

(15) /ir„+i = i^Tn + 7^ (n(/iT„) - /^tJ- 

For the measures /iT„ defined by this procedure, we obtain (using discrete rigorous analogues 
of informal arguments of the previous paragraph) some estimates on the speed with which 
their free energy decreases. This allows us to estimate distances from these measures to the 
set of translates of poo (because they are the only minima of J-"). 

Now, we are taking the true random trajectory p^, and estimate the distance from the cen- 
tered measures p^ to the equilibrium point. To do this at some moment t, we choose an earlier 
moment t', replace the measure fit' by a close smooth measure pc, and consider deterministic 
discrete iterates by flTSl) . On one hand, for this new trajectory the free energy is defined (as 
we have chosen a smooth approximation), so we control the decrease of energy and hence the 
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distance from the centered measure /i^ to poo- On the other hand, an accurate computation 
allows us to control the distance between the random measure /i^ and the approximating de- 
terministic image fit of its smooth perturbation. The sum of these distances then estimates 
the distance from /i^ to poo, and the obtained estimate tends to as t — )■ cx). This concludes 
the proof of Theorem [31 

Finally, to prove Theorem HI one first computes the speed of drift of the center q, and then 
shows that the series of general term |ct„+i — ct„\ converges, and the oscillations osc[T„^T„+i]Ct 
tend to zero. This implies the existence of the limit oi Ct as t ^ oo. 

1.4. Outline of the paper. At the beginning of Section[2], we show the existence and unique- 
ness of solutions to ([2]) starting at any positive moment r > 0. The discussion of this topic 
at t = is postponed to the appendix. In the rest of Section [2], we present some crucial 
preliminary computations which are at the basis of our proofs. Most of the material there 
is not new, except for the combination of stochastic approximation of the empirical measure 
(see [2]) with free energy functionals (see [S]) and the achieving of a bound on the convergence 
rate. Finally, Section [3] is devoted to the proofs of our main results. 

1.5. Acknowledgments. The authors are very grateful to two anonymous referees for their 
useful comments which led to a rewritting of the paper for a better understanding. 

2. Preliminaries 

As usual, we denote by A^(M°') the space of signed (bounded) Borel measures on M'^ and by 
P(R'^) its subspace of probability measures. We will need the following measure space: 

(16) MiR';P):={fieMiR'y, [ P{\y\) \fi\idy) < 00} , 

where |/i| is the variation of fi (that is |yu| := /i"^ + p^ with {fi^,fi^) the Hahn- Jordan de- 
composition of /i: fi = fi'^ — fi~). Belonging to this space will enable us to always check the 
integrability of P (and therefore of W and its derivatives thanks to the domination condi- 
tion ([6])) with respect to the (random) measures to be considered. We endow this space with 
the dual weighted supremum norm (or dual P-norm) defined for p G A^(]R^; P) by 



(17) ||a*||p ■= sup 

>p&C(R'');\ip\<P 



ipdft 



P{\y\)\fi\{dy). 



We recall that -P(|a;|) > 1, so that ||/i||p > |/i(M^)|. This norm naturally arises in the approach 
to ergodic results for time-continuous Markov processes of Meyn & Tweedie [11] . It also makes 
M-iW", P) a Banach space. 

Next, we consider V(R'^;P) = M(R'^;P) n V(R'^). We remark that both M{M.'^;P) and 
V{M.'^; P) contain any probability measure with an exponential tail and, in particular, any 
compactly supported measure. For any k > 0, we also define 



(18) P«(M'*;P):={/iGP(M^P); ||p||p= / P(|x|)/i(da;) < «:}. 

2.1. Existence and uniqueness of solutions. 
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2.1.1. Markovian form; local existence and uniqueness. First step in studying the trajectories 
of ([2]) is to pass to the couple {Xt, Ht). A standard remark is that the behaviour of this couple 
is infinite-dimensional Markovian (and in general, except for W being polynomial, cannot 
be reduced to a finite-dimensional Markov process). This reduction is easily implied by the 
identity 

(19) Ht+s = ^^f^t + —— / Sx^ du. 

t + S t + S Jt 

Note that the second term in the right-hand side of ( fT9l) can be written as j^fiit,t+s], where 
fi[t,t+s] is the empirical measure during the time interval [t,t + s]: 

1 f^'' 

^^[tut2] ■= -. — r / ^^- ^"• 

t2 — tl Jt^ 

Now, passing fit to the left-hand side of flT9|) . dividing by s and passing to the limit as 
s — 7- 0, we obtain the following SDE for the couple {Xt, fit)'- 

.20) / dXi = V2dBt-VW* fitiXt) dt, 

1 At = j{-fit + Sxt). 

For any to > 0, the local existence and uniqueness of solutions to (l20l) . in a neighbourhood 
of to, is implied by well-known arguments: see Theorem 11.2 of J14j . 

However, in order to study the asymptotic behaviour of solutions to ([2]), we should first 
show the global existence of these solutions, in other words, that they do not explode in a 
finite time. It will be done in §2.1.2[ 

Note also that the equation fl2U]) clearly has a singularity at t = 0. To avoid this singularity, 
sometimes the equation fl2U]) is considered with an initial condition {Xr,fir) at some positive 
time r > (and thus for t G [r, cxd)). After the time-shift s = t — r, the system fj2U|) transforms 
to 

.21) / dX, = V2dB, - VW * fis{X,) ds, 

^ ' \ f^s = ^,{-f^s + SxJ. 

In fact, we can restrict our consideration to such situations only (as, anyway, we are in- 
terested in the asymptotic behaviour of solutions at infinity), but it is interesting to show 
that the equation ([2]) has indeed existence and uniqueness of solutions for any initial value 
problem Xq = Xq. It is done in the appendix. 

2.1.2. Center-drift estimates. A natural "reference point" that one can associate to a mea- 
sure fi is the equilibrium point c^ = c{fi) of the potential it generates with W, defined by 
the equation VW * fJ^ic^j) = (see Definition [H §1.2p . that we refer to as the center of the 
measure fi. Also, it will be convenient to consider the centered measure fi'^, obtained from fi 
by the translation that shifts the center to the origin. 

Note that the implicit function theorem allows to estimate (on an interval of existence of 
solution {Xt, fit) to (!20|) ) the derivative q of q := c^^. In particular, we will see that q is a 
C^-function on this interval. 

Indeed, the function (x,t) i— )■ VW * fit{x) is C^-smooth: 

d{VW*fit)ix) = V^W * fit{x) dx + WW * fltix) dt 

= V^W * fit{x) dx + -\/W * {-fit + 5xt){x) dt. 
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and for any {x,t) we have V^PF * fJ't{x) > Cy/I > 0. The imphcit function theorem thus 
imphes that q is a C^-function of t (on the interval of existence of solution), and that 

= j{V'W*fit{ct)y'vW{Xt-ct). 

This implies that the projection of the center drift velocity on the line from q to Xf is 
directed towards Xt, as VW{Xt — Ct) is positive, proportional to Xt — q and 

{{y^w * ^it{ct)Y^ vw{Xt - q), X, - q) > 0. 

This also immediately gives an upper bound on the drift speed: 
(22) |,|<i.£(^^. 

2.1.3. Law of X -center distances: Ornstein-Uhlenbeck estimate. To continue our study, first 
we would like to obtain an estimate on the behaviour of the distance \Xt — Ct\. Namely, we are 
going to compare it with the (absolute value of) Ornstein-Uhlenbeck process, and to obtain 
exponential-decrease bounds on its occupation measure in §2.2.11 

Proposition 1. The process (Xt) can be considered as the first element of the pair {Xt, Zt) 
of processes such that 

i) \Xt-Ct\ <2 + Zt, 
ii) Zt is the absolute value of a 3d- dimensional Ornstein-Uhlenbeck process. 

Proof. From 

dXt = V2dBt - VW * fitiXt) dt 

Ct = i {V^W * fitict))-' VWiXt - Ct) 

- Ct\, while it is non-zero, satisfies the SDE 

,dBt) + .^'^ . dt 
J \Xt-Ct\ 



one 


obtains 


that the difference 


Xt- 




d|X,. 


-ct| 


^ V2( 
(Xt 


Xt- 
Xt- 

-Ct 


-Ct 
-Ct 

- VI 



, viy * iit{Xt) + ^(v'ly * iit{ct)Y^^w{Xt - ct)\ dt. 



\\Xt-Ct 
In the same way, the desired Zt shall satisfy the equation 

(23) dZt = V2d^t - {^Zt - ^) dt, 

where 7 is also a Brownian motion. So, take a one-dimensional standard Brownian motion /5 
independent of the Brownian motion B and let 7 be defined as 

(24) d7t = a{\Xt -ct\) f ^^^, dB^ + ^l - a^Xt- ct\)d[5t, 

W^t — Ct\ J 

where a : [0, -|-oo) — )■ [0, 1] is a C°°-function which is identically zero in some neighbourhood 
of and a{r) = 1 for any r > 1. The process Z is then defined by (|23|) . 

We point out that, as B and /3 are independent, i? is a d-dimensional Brownian motion while 
(3 is 1-dimensional. It follows that Z defined by fl23|) is the absolute value of a 3(i-dimensional 
Ornstein-Uhlenbeck process. 
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On the other hand, for any t, either \Xt — Ct| < 1 and then automatically \Xt — Ct\ <2 + Zt, 
or |Xt — qI > 1 and then both |Xi — Ct| and Zt share exactly the same Brownian component 
(as a = 1), with the inequality between the drift terms oi 2 + Zt and \Xt — Ct\: 

(25) -^Zt + ^^>-Cw\Xt-Ct\ + j^—,> 

d 



V \Xt-ct\J 



d~ 


-1 


\Xt- 


-ct| 


Xf 


-Ct 



\Xt-Ct\ 



^{y^w*^^t{ct))-^vw{Xt-ct] ^* ""* 



t \Xt-Ct\ 

as soon as \Xt — Ct\ > -^j^Zt- A comparison theorem concludes the proof. D 

2.1.4. Global existence. 

Proposition 2 (global existence). For any r > and for any initial condition {Xr,Hr), the 
solution to (120]) exists (and is unique) on the whole interval [r, +oo) . 

Proof. As we already have the local existence and uniqueness, it suffices to check that the 
solution Xt cannot explode in a finite time (this impossibility will imply that the measures fit, 
as the occupation measures of Xt, also stay in a compact domain for any bounded interval of 
time). 

Let us introduce the increasing sequence of stopping times tq = and 

Tn := inf {t > r„_i : \Xt\ > n} . 

In order to show that the solution never explodes, we use the comparison of Xt — Ct with the 
Ornstein-Uhlenbeck process Zt (see §2.1.3p . So, we have for the corresponding Z, that 

|^min(t,T„) ^ Cjnin(t,T„)| S ^ + ■^min{t,r„)- 

As Z does not explode in a finite time, letting n go to infinity, we conclude that Xt — Ct does 
not explode in a finite time. To conclude, one has to use the inequality f l22|) : 

t G]y t G]Y t G\Y 

Any trajectory of Z being bounded on any finite interval of time, the integral J ds is 

finite for any t > r. So, the process {Xt,t > 0) does not explode in a finite time and there 
exists a global strong solution. D 

2.2. Exponential tails estimates. 

2.2.1. Estimates for the centered empirical measure. We shall now estimate the behaviour of 
the centered measures /i^. Namely, we are going to prove that these measures are exponentially 
decreasing. For shortness and simplicity, we introduce the following 

Definition 2. Let a, C > be given. Then 

(26) Klc ■■= {/iG^(K'); Vr, /i({y; |y| > r}) < Ce-"^'}, 

(27) K^,c := {/xG^(K'); f^"" ^ Kc}- 

Also, for non-probability positive definite measures, we denote the spaces defined by the same 
inequalities by K'^q and Ka,c- 



ERGODICITY OF SELF-ATTRACTING MOTION 13 

For what follows, we need one easy lemma. 

Lemma 1. Let Z be the absolute value of a 3d- dimensional Ornstein-Uhlenbeck process. Then, 
there exists Ci > 0, such that for almost any trajectory Zt, one has almost surely 

3T : Vt > T, Vr > -\{s<t: Zs> r}\ < Cie-\ 

Proof. Note that the Ornstein-Uhlenbeck process is ergodic, with stationary measure 70(7 = 
g-Cwkl /2_ fjj-^g function /(x) = e'^' is 7oi7-integrable. Hence, by (Birkhoff) ergodic theorem, 
almost surely 



1 '•* 



fiZ^s ^ / f{x)d^ou{x) =■■ I. 



Thus for all t large enough, |/q e'^^'ds < / + 1. Applying Chebychev's inequality, we see that 
for all r > 0, 

^|{s<t:Z, >r}|<(/ + l)e-^ D 

The main result of this subsection is the following, showing that the measure /ij belongs to 
the set Ka^c- 

Proposition 3. There exist two constants a,C > such that a.s. at any sufficiently large 
time t, Ht e Ka^c- 

To prove this proposition, we need two intermediate lemmas, which proofs are postponed. 

Lemma 2. There exist ao,Co > such that a.s. for any sufficiently large time t, ^[t/2,t]{' + 

ct/2)e^S„,co- 

Lemma 3. Let a^, Cq > Q be fixed. Then there exist a, C, C such that the following holds. 
Assume that there are given a measure /i G V{M.'^\P), a measure u{- + c^) G -K^qq Cg o-nd a 
coefficients < A < 1/2. Then if 11 can be decomposed as ft = fi^^^+fi^'^^ with fi^'^\- + c^) G K^c 
then for the decomposition (1 — A)/i + Az/(- + c^) = (1 — X)fi^^^ + ((1 — A)/i*^^^ + Xh'{- + c^)), 
one also has ((1 - X)fi^'^^ + Az/(- + c^)) G K^^c- 

In other words, this lemma provides an "induction step" for showing that "a big part of 
the centered measure has exponentially small tails" for a procedure of repetitive mixing with 
measure with exponential tails (this not being obvious, as the center can be shifted by such a 
procedure). 

Proof of Proposition^^ First, let us estimate the drift of the center. Namely, taking to- 
gether (1221) and Proposition [H we have 

\ct\ < T^P{\Xt - qI) < -^P(2 + Z,) < ^P{Z,), 

l(^W i^'^W i^'^W 

for the corresponding Ornstein-Uhlenbeck trajectory Zt. 

On the other hand, Z is a Harris recurrent process and P{Z) is integrable with respect to 
the Gaussian measure, thus due to the limit-quotient (or Birkhoff) theorem, almost surely 
there exists a limit 

hm I [ P{Z,)ds = [ P{\z\)d-fou{z) =: /. 
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So, almost surely from some time ti we have 

Vt >ti, - / P{Zs)ds<I + l. 
t Jo 

Therefore, after this time we can estimate the displacement of the center between the moments 

t/2 and t: Wt > ti 

\ct/2 - qI < / |c.| ds < / -P{Z,) ds<^ [ P{Z,) ds < 2C{I + 1) =: C3. 

Jt/2 Jt/2 S t/^ Jo 

In fact, the same estimate holds for any t' between t/2 and t: 

\ct' -ct\ < C3. 
This immediately implies that for any t > ti and n G N such that 2~"+^t > ti, one has 

\Ct - Ct/2"| < C-iU. 

Let us now apply Lemma [3l First let us decompose, for any t G [ti,2ti], the measure 

/i2t as \iit + i/^[t,2t], then the measure /i4t as i/ii + (|/i[t,2t] + |yU[2t,4t]), • • •, and finally the 

measure /i2"t as -^iit + {■^l^[t,2t] + ■ ■ ■ + |yW[2"-if,2"i])- An induction argument, together with 

Lemma [2l immediately shows that in each such decomposition, the second term shifted by the 

corresponding c{ix2h) belongs to K^^j. The only part that is left to handle is ^/xj. But the 

distance between q and C2>h does not exceed C^n, and the centered measure nl is compactly 

supported. So it is contained in a ball of some (random) radius R that can be chosen uniform 

over t G (ti, 2ti). Now the measure -^^t is of total weight 2~" and it vanishes outside a radius 

R ball. If a is small enough so that e"*"^ < 2, then for any r > C^n + i?, we have 

1 1 

l^^^t{\y - C2"t| > r) < — /ii(|l/| > r - C^n) = 0, 

and for r < C^n + R and n big enough, 

^f^t{\y - ca-tl > r) < 2-" < e-""^3e-ai? < g-ar_ 

The middle inequality comes, for n large enough, from a comparison between exponent bases, 
goCs ^ 2^ with respect to which a multiplication constant e"""^ is minor. Finally, joining the 

obtained ^/it(- + C2^t) e -R'° ^ and (^/i[4,2t] H h |/X[2"-ii,2"i]) (■ + C2"i) G -ft'^ ,^, we obtain 

Proof of Lemma [H This lemma immediately follows from Lemma [H once we notice that 

2 

/^[t/2,i](|y-ct/2| > ^) = - |{s : t/2 < s < t, |X^ -Q/al > r}| 

2 I I 

< - |{s : t/2 < s < t, |X, - c,| > r - |q/2 - c,|}| 

< -|{s: s<t, Z, >r-C3}| <Coe"''^^-e-°o". D 

Proof of LemmalBi First, let us estimate the position of the center of /i in a way that is linear 
in A and does not depend on a and C — thus in particular, proving the statement i). Indeed, 
cp_ is the minimum of the function W * fi. At the point c^, the gradient of this function can 
be bounded as 

I VW^ * j2\c,\ = 1(1 - X)VW * /x|,^ + WW * u^ < XM- + c^)\\p < C'X, 
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because the norm ||z/(- + c^)||p is uniformly bounded due to the condition z/(- + c^) G Ka^^^Co- 
Now, restricting the function W * jl on the hne joining c^ and c^, that is considering 



f[s) = W *^i[Cf, + s 



\r~ — r 
I ^M M I 



one sees that |/'(0)| < C'X, /'(|c^ — c^|) = 0, /" > Cw, what imphes |c^ — c^| < ^A. 

Let us now estimate the measure ((1 — A)//*^^^ + Az/) (|y — c^| > r). Indeed, note that 
{y ■ \y- C/i| > r} C {y : \y-c^\>r- C"\}. Thus 

/i(|l/-c(/i)|>r) < f,(\y-c{ix)\>r-C"\) 

< (1 - A)M(2)(|y - c(^)| > r - C"A) + An(^^)(|i/| > r - C"A) 

< (1 - A)Ce-"('^-^"^) + ACoe-"°(^-^"^) 

(28) < e^(^""-^/2)Ce-°" - A f-e^»'^-"> - Coe"°^"M 6""°^ 

Once a is small enough so that C"a < 1/2, a < a^ and once C is greater than 2Coe"'''" , the 
right-hand side of fl28|) is not greater than Ce""*", what concludes the proof. D 



2.2.2. Estimates for the centered measure U. 

Lemma 4. For any n > 1, the map U restricted to Vk{^'^] P) is bounded and Lipschitz. 

Proof. First, we need to show that Z{fi) is bounded from below on Vk{^'^', P)- For /i G 
Pk(M'^;P), the domination condition (EJ implies that W * /i(x) < ||/i||pP(|x|) < kP{\x\). So 
we have: 

Z(/i) = / e-^^'^^^Mx > / e-"^(l^'l)dx. 



Now, using that W*fi{x) = J W{x-y)fi{dy) > ^ / \x-y\'^fi{dy) > ^ / (^ - |yp) ^{dy) > 
~2' ( 4 /t j , we hence have the following bound for 11 (/i): 

(29) ||n(/x)||p< f/ e-"^(l"l)dx) ■ / PdxDe-^^l'I'/^-'Ma; =: a. 

Note that 11 is C^ on 7'(]R'^; P) endowed with the strong topology. As the set of probability 
measures has no interior point, we have to specify the meaning of C^: there exists a continuous 
linear operator DIi{^) : AioiM.'^^P) — )■ Aio(W^; P), continuously depending on /i, such that 
\\U{fi')-U{fi)-DU{fi){fi-fi')\\p = 0(||/x-/x'||p) provided that /x' G P(R'^; P) and /x' converges 
toward fi. Indeed, it is easy to see that 

Dn(/i)-z/ := -(W^*z.)n(/i)-^^|^4^e-^*^ 



Z(^) 



-{W * u)U{fi) + I W*u{y) dy 



Z(/i) ^ Z(/i) 
(30) = -(w*iy- f W*iy{y)n{ix){dy)jU{fj,). 
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Now, note that the norms ||-Dn|| are uniformly bounded for fi G Vk{^'^] P) (for any given k). 
Indeed, fix z/ G Mq{W^;P). Since \W * v{x)\ < ||z/||pP(|x|), we find that 

||Dn(/i)-z/||p<(i + a)||//||p / P'(|x|)n(/i)(dx). 

For yU G Vn{W^\P), the same computation used for the bound (1291) on the norm of n(/u) 
enables to control the last integral. Hence, we deduce a bound (call it C^) on the norm of the 
differential. Thus, 11 is Lipschitz as stated. D 

We prove now the exponential decrease for the centered measure n(/i). 

Proposition 4. There exists a positive constant Cn such that for all fi G V(R;P), for all 
R>0, we have n(/i)(|x - c^\ > R) < Cue'^"^^. 

Proof. Note first that, imposing a condition Cn > e^*"^', we can restrict ourselves only on 
R > 2: ioT R < 2, the estimate is obvious. The measure n(/i) has the density -^^e~^*^^^\ 

To avoid working with the normalization constant Z{fi), we will prove a stronger inequality, 
that is 

n(/i)(|x - c^\ >R)< Cne-^^« ■ n(/i)(|x - c^\ < 2), 
which is equivalent to 

Pass to the polar coordinates, centered at the center c^: we want to prove that 

/ / e-^*'^(^^+^'')A'^-MAdt; < Cne"^'^^ / / e-^^'^^^^+^^U'^-MAdt;. 
Js'i-^ Jr Js''-^ Jo 

It suffices to prove such an inequality "directionwise" : for all v G E>'^^^, for all i? > 2 

oo ^2 






But from the uniform convexity of W and the definition of the center, the function /(A) = 
W * ji{c^ + Xv) satisfies /'(O) = and Vr > 0, f"{r) > Cw Hence, / is monotone increasing 
on [0, oo), and in particular. 



(31) / e-^^^U-^-MA > e-^(2) / A'^-MA =: C^e-f^^\ 

Jo Jo 

On the other hand, for all A > 2, /'(A) > /'(2) > 2Cw, and thus /(A) > 2Cwi\ - 2) + /(2). 

Hence, 

(32) 



'R JR 

Comparing fl3T]) and fl32|) . we obtain the desired exponential decrease. D 
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2.3. A new transport metric: Tp-metric. Usually, to estimate the distance between two 
probability measures, one introduces the Wasserstein distance. Indeed, for /ii,/U2 G V(M.^; P), 
we define 

iy|(/ii,/X2):=inf{E(|ei-6l')}, 
where the infimum is taken over the random variables such that {law of ^1} = fii and {law 
of (,2} = /^2- In our setting, for a measure /x, the corresponding probability measure n(yu) is 
defined using the convolution W * fi. So, it would be rather natural to use a distance, looking 
like the one for the weak* topology, but allowing to control W * fi for our unbounded function 
W. This motivates to introduce a new metric looking like the Wasserstein distance: 

Definition 3. For /ii,/i2 € V{M.'^; P), we define the P -translation distance between them as 
(33) rp(/ii,/i2):=inf|yy" P{\f{s,co)\M{s,u)\dsdF^ 

where the infimum is taken over the maps / : [0, 1] x fi — )• M, where f2 is a probability space, 
such that {law of /(O, ■)} = /ii, and {law of /(I, •)} = /i2. 

We also denote the 7p-distance between two c-centered measures by T^(/ii, /i2) = 7p(/ii(- + 
c),/i2(- + c)). 

Remark 5. In dimension one, we have the equivalent definition: 

Tp(/il,yU2) := / P{\x\)\l2i{{-00,X]) - IJ,2{{~00,x])\dx. 

Jm. 
The following lemma will be useful to show the convergence of the empirical measure in the 
Vr2-nieaning, as Proposition E] shows. 

Lemma 5. Let //i,/X2 G V{M.'^] P). There exists a constant C > 0, such that 

1^2^(yUi,/i2) < CTp(/il,/i2)- 

If moreover ^i, ^2 belong to a set Ka^Co, ^^^'^ there exists C > 0, such that 

rp(/ii,/i2)<c"iy|(/^i,/^2). 

Proof. Suppose that /Ui,/i2 G Ka^Co- Take ^1,^2 realizing the optimal VF2-transport between 
them, and let us estimate the Tp-cost of the same transport. Indeed, 

7p(/ii,/i2) < / |6-6l^(max(|^i|,|^2|))dFi„ax(6,6) < W^2^(/^i,/^2) / P^(max(|^i|,|^2|))dFmax(6,?2)' 
where the second inequality is the Cauchy one. As /ii,/i2 G Ka,Co, ^^ conclude that 



/"p^(max(|ei|, |6l))dF^ax(a,6) < J P\r )dmax{0, 1 - 20^6- 



--: a < cx). 



Let now ^1,^2 be two random variables corresponding to the 7p-optimal transport of /ii to 
yU2- We then have 

W^|(/^i,/^2) = / 1^1 - 6pdP < / 1^1 - 61 ■ 2max(ei,6) < / 16 - 61 ■ ^^"'^^^f ^^^^^^ 

(34) <Crp(/ii,/i2). 

Indeed, the inequality (IMl) is due to the fact that the path between ^1 and ^2 either stays 
outside the radius max(|,^i|, |^2|)/2-ball centered in 0, in which case we estimate its length 
from below as |^i — 61; or this path has a part joining the maximum norm vector to this ball, 
which is of length greater than max(|^i|, |6l)/2 > |6 ~ 61/4- □ 
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It is clear from the definition that 7p is a distance; and also taking into account that 
\P'\ < P, one easily has 

(35) ||/^2||p < ||/Ui||p + Tp(Aii,/i2). 

Thus, the set V{R'^] P) is 7p-complete. Indeed, a Tp-Cauchy sequence (/i„) will have a weak 
limit fi and it is easy to check that ||/i||p = lim ||/in||p < oo. So, ft G P(M°';P). Now, we 

are going to estimate the deviance of trajectories in terms of 7p-metric, a result that will be 
useful in §3.1[ 

Lemma 6. For yUi,/i2 G V{M.'^; P), the following statements hold: 

1) The map c is locally Lipschitz in the sense of Tp-metric: 

|c(/ii) -c(/i2)| < -^mm{P{\c{fii)\),P{\c{fi2)\)) ■ Tp{fii, ^12)] 

2) For all V G R'^, we have Tpifi, fii- + v)) < \v\P{\v\)\\fi\\p; 

3) There exists Cp > such that 

^c, ^^ P(x+H) , /(l + Cp|t;|)-rp(/i,z/), |t;|<l 

r^/i,i/< sup — Tp{fi,u) < ^ 

4) For all k>0, fi^ : V^{M'^; P) -^ V{R'^; P) is Tp-Lipschitz. 
Proof. 1) Denoting by Ci (resp. C2) the center of /ii (resp. /i2), we have 

VW * /i2(ci) = VW * /ii(ci) + VW * (/i2 - /Ui)(ci), 

thus |Viy*/i2(ci)| < P(|ci|)||;U2 — /ii||p. Joining the points Ci and C2 by a line, recalling that 
due to the uniform convexity of W, the second derivative oi W * fi2 along this line is at least 
Cwf^2{R^) and noticing that VW * /X2(c2) = 0, we obtain 

(30) |C2-Ci|<— ||/i2 -/Ul||p- 

Similarly, we have |c2— ci| < ^'^^ ||/^2— /^i||p- So, the result follows as ||/i2— yUi||p < 7p(yUi,/i2). 

2) We have by definition of 7p that 

rp(/i,M- + l^l))= / Kdx) [ Pi\y\)dy< [ \v\P{\x\ + \v\)fiidx)<\v\Pi\v\) [ P(|x|)/i(dx) 

3) For any transport /(s,^) between /i = {law of /(0,Ci;)} and u = {law of f{l,u)}, the 
map f{s,uj) — V is a transport between fi'^ and u'^ of price 



JnJo x>o P{x) JnJo 



Pi\fis,u)\)\f:is,u)\dsdPiu). 



The left-hand side is an upper bound for T^(/i, i') and passing in the right-hand side to the 

P{x+\v 
c>0 p(x) 



infimum over all the possible transports /, we obtain the desired sup^>Q pj^v' 7p(/x, u). 
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4) Suppose that fi,!^ E V^iM'^', P)- Then, by the preceding points, we have 

Tp(yU(- + C^),Z/(- + C^)) < Tp(yU(- + C^),Z/(- + C^)) +Tp(z/(- + C^),Z/(- + Ci,)) 

< Pilc^DTpifi, u) + \c^ - c,\Pi\c^ - c,|)||z/(- + C,)\\p 

< Pi\c,\)rpifi,iy) 

+ P(|c^-c.|)^min(P(|c^|),P(|c,|))||z/||prp(^,z/). 

Remark that, as /i, z/ G Vk.(J^'^', P), the norms \c^\ and |c^| are uniformly bounded, as well as 
||z/||p, thus 

TpiK- + c,), K- + c.)) < (P(|c,|) + Pile, - c.|)^ min(P(|c^|), P(|c.|))||z/||p^ TpifJ., u), 

where P(|c^|) + P(|c^ — Cy|)^min(P(|c^|), P(|c,^|))||z/||p is uniformly bounded by some con- 
stant Ck, which is the Lipschitz constant. D 

2.4. Free energy functional. We recall from §1.3.41 that the free energy of a measure is 
defined as 

^■(/i) = nifi) + 2 // K^Wix - y)fi{y) dxdy, n{fi) = / fi{x) \ogfi{x)dx. 

The free energy of a non-self-interacting gas in an exterior potential V is defined as 

•^y(/^) = '^(a*) + / iJ.{x)V{x)dx 

and the map 11 associates to a measure fi the probability measure ^e~^*^^^^dx (when W * fj, 
is well-defined). 

The first auxiliary statement implies that, as we mentioned it in §1.3.41 11 (/u) is the unique 
global minimum of J-'w*ij.- 

Lemma 7. For any potential V such that e~^ is integrable, the probability measure Z~^e~^ 
is the unique global minimum oj Ty on P(R°'). 

Proof. Let n = Z^^e^^ . Then, for any arbitrary absolutely continuous measure p, letting 
p{x) = Ze^^^^u(x) be its density with respect to /x, we see that 

JV(i^) = / (V"(x)+logz/(x))z/(dx) = / (logp(x)— logZ)z/(dx) = / p(x) logp(x)/i(dx)— logZ, 

jRd JR'i jRd 

and thus Jensen's inequality, for the convex function plogp, leads immediately to the conclu- 
sion, n 

Now, for the free energy functional, McCann [T0| proved the following 

Proposition 5 (McCann). There exists a centered symmetric density p^, which is a unique, 
up to translation, global minimum of T . Moreover, T is a displacement convex functional, 
that is for two probability measures po, pi and the W as ser stein- optimal transport between them 

is = {l- 5)^0 + s^i, 

where po = {law of C,q}, pi = {law o/^i}, E|^o ^ GP = W^2^(/io,Pi); one has 

J^{{law o/^s}) > (1 - s)J^{po) + sJ'ipi). 
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Finally, the transport distance from a centered measure fi to poo can be estimated as 

(37) W^|(/i,Poo)<7?--F(/i|Poo), 

where J'(p|poo) = -7^(p) - ^(Poo)- 

Remark 6. i) The uniqueness of the minimum comes from the strict displacement convexity 
of the restriction to the space of centered measures. 

ii) The functional J-' is not convex in the usual sense, due to the self-interacting part. 

iii) Inequality flT^ together with Lemma [7| immediately imply that the minimum of J-" is 
also a fixed point of II. 

Finally, as we are going to work in §3.1.3l with the discretized flow, we will need two auxiliary 
statements for the free energy: 

Lemma 8. For all absolutely continuous measures /x, z^ G V{M.'^]P) of finite free energy and 
for all X G [0, 1], we have 

(38) J'((l - A)p + Az/|poo) < ^(pIPoo) - K^M - 'f^^{^)) + 

+ Y / / (/" ~ u){x)W{x - y)(p - u){y) dxdy, 

where (Pfj.{-) := J^w*fi{-) is the free energy in the fi-generated potential. 
Moreover, for all absolutely continuous p G V(M.'^; P), we have 

(39) v.^(p) - ^^{v) = ^(p) - ^(z/) + ^ jjifi - v){x)W{x -y){fx- u)iy) dxdy. 

Proof Note that H{{1 - A)p + Xu) < (1 - A)H(p) + AH(z/) = H(p) - A (H(p) - n{u)). So, it 
suffices to prove (1551) with entropy terms removed form both sides (from both J-" and y^^ in the 
right-hand side). After this removing, the formula becomes a Taylor expansion for a degree 
two polynomial. The same holds for (!39|) . with a remark that the entropy terms are exactly 
the same in both sides. D 

3. Proofs 

3.1. Proof of Theorem [3l In fact, we will prove a stronger statement, controlling the speed 
of convergence in the sense of the transport distance: 

Proposition 6. There exists a > such that almost surely, as t ^ oo, 

rp(p^,Poo) = o(e-'^'="^, 
where k is the degree of the polynomial P, as well as 

W^2(/i?,Poo)=o(e-«'="^ 

The proof of this statement will be decomposed into several propositions. We first present 
them all, postponing their proofs; then deduce from them Proposition [61 Finally, we prove 
these propositions. 

In order to prove this statement, as it was announced in §1.3[ we will discretize the random 
process. Namely, we define the sequence T„ of moments of time as T„ := n^/^ and then. 
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I/O 

AT„ := Tn+i — Tn is of order TJ ■ Also, for what follows, we will associate to a random 
trajectory {Xt,t > 0) thz sequence (L„) defined as 

(40) L„ := max \Xt - CtJ < C^logT„. 

0<t<T„+i 

An easy conclusion from the Ornstein-Uhlenbeck comparison §2.1.31 and logarithmic drift of 
the center is that almost surely L„ < C3 log n and L'„ < log n for any n large enough. 

Now, let us state the first of the propositions mentioned above, the one allowing to estimate 
the "Euler-method" one-step error in the description of the behaviour of measures /it: 

Proposition 7. Almost surely there exists uq such that for any n > uq, we have 

where (3 = min (SCw, ■^) ■ 

Associated to the moments of time T„, consider the following, roughly speaking, Euler- 
approximation maps for the flow rh = j(Il{m)—m), with the knots chosen at the moments T„: 

Definition 4. For any i < j, define ^ : P(M°', P) -^ V{R'^, P) as 

AT- 

^ = zd, $r(M)=/^ + 7^(n(/x)-/^), $H$^-io---o$p. 

Let us first exhibit an invariant set for $. 

Lemma 9. For any a,C as in Lemma 0, corresponding to oq = Cw a'^d Cq = Cu (from 
Proposition^^ , if ft E Ka^c and i < j, then $i(/i) G Ka,c- 

Proof. This is a direct corollary of Lemma |3l D 

Denote, for a probability measure ft and for a number h > 0, hy jj,^^'^ the "smoothening 
convolution" 

where Uh{0) is the radius h ball in Mf^, centered at the origin. 

The following proposition allows to compare the deterministic Euler-like behaviour of the 
smoothened, at some moment Tj, measure with the true random trajectory: 

Proposition 8. There exist some constants A, 6*1,6*2, C3 > such that almost surely there 
exists no for which the following statements hold. For any j > i > uq and any h > 0, 

(41) rr^($K4?),/^T,) < i:^(Ar.)-^ (i)^ + ^^' (1)^' 

k=i ^ ^ 

provided that the right-hand side of ( BTj) does not exceed 6*3. ^4/30, under the same condition, 

Next, we have to show that the deterministic trajectory of an absolutely continuous measure 
sufficiently fast approaches the set of translates of poo- To do this, due to the estimate ( 1H7|) . 
it suffices to estimate the free energy: 

Proposition 9. Let fi G K^^. Then, there exist ai,C*4,C*5 > such that almost surely there 
exists Uq for which the following statements hold for any j > i > uq: 
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i) ^($K/")IPoo) < C4 + |(^(/i|poo) - C4), 
ii) ^(<l'K/i)|Poo) < Cse-"^ ''-V^o.mn) ,f:p^^\p^) < 2C4. 

Now, modulo these propositions, we are ready to prove Proposition [6l 

Proof of Proposition\^ Recall from Proposition [7] that /3 = min(8Ciy, (5(i)~^). Note first that 
the distances T^"^" {fit, fiT„) for t G [Tn,Tn+i] are uniformly bounded for n sufficiently big by 






where L„ is defined by fj40|) . Hence, it suffices to check the estimate for the subsequence of 
moments T„: 

Now, for any sufficiently big n, take i := [n^"*^], where a small 5 > will be chosen and fixed 
(in a way that does not depend on n) later. Then, considering for some h > a. smoothened 
convolution /x^, and its Euler-image $"(/Xj, ), we have by Proposition [8] 

(42) 7^-(PT„,$r(/^ff)) < E^^^ (p^)%c7,/. (|)^ 

provided that the right-hand side does not exceed C3. 

Denote by const a generic constant. Let us estimate the first term in the right-hand side: 

n—l . rn /m \ A n—1 



AT /T \ 9 /T \ 

(43) E|^(A^o-(|) si:|(An)-(|) ^ 

E-(Ar.)-^(jf) < const .„L_L^_j 



(l + A+|)<5-f 



So, for any fixed choice of 5 < j-^^tow:, the first term in the right-hand side of ( 142|) will 
decrease as a negative power of n and thus quicker than e~° viogTW_ 

Take now /i = ^ ( ^ j . For such a choice of h, the second term in the right-hand side 
of fH2D is not greater than ^ ~ n^^. So it also decreases quicker than e~" viogT^ and thus 

Finally, we have to estimate 7^^"($"(/i^. ),poo(- + ct„))- To do this, it suffices to estimate 
the free energy ^^(^^(/ij-. )), as 

Tp ((<l>^(pi.t^))^Poo) < const W^l ((<|.^(pi.t)))^Poo) < const ^(<l>r(pi.^))). 
Indeed, 

J-(Aig)) = ?^(/.g)) + fL^f^{dx)W{x - y)/xgHd|/). 



The first term here does not exceed —\ogvol{Uh{^)) < d ■\ \og{h/d)\ (as the density of p^.. 
does not exceed {h/d)~^), while the second term is bounded. Thus J-'(pj.. ) < Celogn for 
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some constant Cq. Hence, from the first part of Proposition [9l for j = 



r 1 V^^ ' 
^ log n 1 I 



'J 



sj7„W 



Applying the second part, with $:^(/i^/) as a starting measure, we obtain 



Thus, J^(<l>^(;u5.^^)) < e-"'^^/i°g^ and hence 

Tp [mif^^fiW, Poo) < const W^ (($^(/^g^))^ Poo) < const e"" ''^^^^^^. D 

Let us now prove Propositions [7H91 

3.1.1. One-step error estimate. This section is devoted to the proof of Proposition [71 

To estimate the difference between the occupation measure of Xt on [T„,T„+i], and the 
measure Il{^t), "we will first introduce another process, for which n(/i2-„) is the stationary 
measure: the process with "frozen" measure /iT„- More precisely, on [T„,T„+i) we consider a 
process Y with some choice of Yt„ , satisfying 

(44) dVt = V2dBt-VW* iiT„ {Yt) dt, 

generated by the same Brownian motion Bt as Xt. In other words, the couple (Xj, Yt) satisfies 

...x f dXt = V2dBt - WW * fitiXt) dt 

^ ' \dYt = V2dBt - VW * firAYt) dt. 

The following lemma allows to control the difference between them: 

Lemma 10. For all t G [T„, Tn+i] we have 

AT 

(46) \Xt - Yt\ < e-^-(*-^")|XT„ - YtJ + --^P(2L„). 

Proof. The process Xt — Yt is of class C^. We compute 

^(Xt - Yt) = -{VW * fitiXt) - VW * firAYt))- 
dt 

Adding and substracting VW * fiXni^t), we see 

d{Xt -Yt) = - [VW * {fit - PrJiXt) - {VW * fiT^Yt) - VW * fiT^Xt))] dt. 
The last term can be rewritten as 

-{VW * /iT„(r,) - VW * f^rAXt)) = ^ f^ f V2iy|.y,+(i_.)x.-x. ■ {Xt - Yt)dud.s. 

J-n Jo Jo 

Noting the first term as Dt, and putting a scalar product with Xt — Yt, we see 
- — iXi-Fip = {Dt,Xt-Yt) 

- yJ,\J, V'^|(„y.+(i-„)x.)-xldu ■ {Xt-Yt),Xt-Y}j ds 
< {Dt,Xt-Yt)-Cw\Xt-Yt\\ 



24 VICTOR KLEPTSYN, ALINE KURTZMANN 

Thus, ^|Xi - Ftp < -2Cw\Xt - Yt\^ + 2|A||^t - YtV Redividing by 2\Xt - Yt\, we obtain 

(47) ^|Xi - Y,\ < lAl - Cw\Xt - Yt\. 

Finally, notice that |-Di| < P{2Ln)^^, as it is the difference between the forces generated at 
Xt by /iT„ and by fit = HTu + ^^ {l^[TnA - l^Tr)- Solving iit = P(2L„)^ - CwUt, we obtain 
the desired estimate for the difference \Xt — Yt\ on the interval [T„,T„+i]. D 

For what follows (see Proposition [TUJ and Lemma [T^ below) . we will have to assume that 
the initial distribution of Yt„ is absolutely continuous with respect to 11 (/i^^), and to use an 
estimate on its density. So finally, we define the process Y^ for all t in the following way: for 
every interval [T„, T„+i) the initial value Yt„ is chosen randomly with respect to the restriction 
of n(/iT„) on the unit ball Wi(ct„). On each new interval, the choice is independent of X and 
of all the past. Then, inside the interval (T„,T„+i), the couple (X„,y„) satisfies fHSj) . 

Let us compare the occupation measures of the processes X and Y on these intervals of 
time. Denote by /i[T„,T„+i] the occupation measure of Y on the interval [T„,T„_|_i]. Then, we 
have the following: 

Lemma 11. For any family of choices Yt„ € Ui{ct„), we have 

7^^" (/^[T„,T„+i], /i[T„,T„+i]) = o{T-^/^), as n -> oo, 
provided that for n sufficiently big Ln < C'^ log ra. 

Proof. The measures fi[T„,T„+i] and /i[r„,r„+i] are both images of the normalized Lebesgue 
measure -^rLeh[T„,Tn+i] under the maps X, and Y, respectively. So, consider the transport 
^s(t) = (1 — s)Xt + sYt between them. 
Using this transport, we have an estimate 

1 fTn + l /•! 

r;""(/i[T„,T„^,],/i[T„,T„^,]) <^jr / ^((1 - s)Xt + sYt - ctJIX, - Yt\ dsdt 

1 pTn + l 

(48) -M^ P(max(|X,-CTj,|F,-CTj))|X,-r,|dt. 

By definition of L„, we have Vt G [T„,T„_(_i], \Xt — ct^\ < L^ and due to Lemma [TOl 

AT 

\Yt-Xt\ < e-^-(*-^")|Xr„-rrJ + --^P(2L0 

AT 
< L„ + 1 + --^P(2L„) < L„ + 2, 

provided that Ln < C3 log n and n is sufficiently big. This implies that 

\Yt - ctJ < \Yt - Xt\ + \Xt - ctJ < 2L„ + 2. 
Now, substituting the obtained estimates to the right-hand side of (j^6l) . we see that 

'7p^"(/^[T„,T„+i],/i[T„,T„+i]) < 

1 /"'^"+i / AT \ 

- A^ X ^^^^" ^ ^^ ■ (^"''"^"'^"^^^" + 1) + ^^^^^"V "^^ - 

^ P(2L^ + 2)P(2L„)AT^ ^ P(2L„ + 2)(L„ + 1) ^ ^.y-l/5^ 
CwTn Cvi/ATji 
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1/3 

we have used that AT„ ~ T„ , and once again the logarithmic growth of L^ 
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D 



Now, we will compare the occupation measure /i[T„,T„+i] with Il{fiT„). To do this, we use 
Proposition 1.2 of Cattiaux & Guillin [6] (see also Wu [17J), stating that the trajectory mean 
of a function tp is, with a probability close to 1 that can be exponentially controlled, close to 
its stationary mean. Namely, this proposition says the following: 

Proposition 10 (Cattiaux & Guillin |6]). Given a process ^ with a stationary measure m and 
Poincare constant Cp, an initial measure v and a function ip satisfying \iIj\ < 1, one has for 
any < p < 1 and t > 



F„|i 



^Pi^s) ds — ^jjdm > p ) < 



du 



dm 



exp 



L2(m) 



tp' 



iCpVarmW 



We will use this proposition with ip being the indicator function tp = 1m of various sets M: 
it then allows to compare the occupation measure of the set M to its n(/iT„)-iiieasure. 

We know that m = Il{fiT„) is the unique stationary measure of the drifted Brownian mo- 
tion fj44|) . Also, the Poincare constant for this process is 2Cw (see [1]). 

To proceed, we have to declare the initial measure i/ = i/„ for Yt„, and we choose it to be 
the measure II{ht„) restricted to the ball Ui{ct„) and then normalized accordingly. Then, 



dz/„ 



dn(/XT„; 



n(/iTj(Wi(cTj) 



< Ce = const. 



i2(n(MT„)) 

the latter inequality is due to the exponential tails of n(/iT',J. Having made these choices, we 
are going to prove the following 

Lemma 12. As n ^ oo, we have almost surely 

'7^""(/i[T„,T„,,],n(/iTj) = o{{AT„y-H^^-'^.)). 

Proof. The previous estimates imply that the process Yt on [T„,T„_|.i] almost surely for all n 
sufficiently big stays inside the ball Up^X^Tn)) where i?„ := 3L„. Now, take this ball and cut it 
into some number Nn parts Mi, . . . , M^^ of diameter less than Sn '■= ^^ (by cubic the grid 

with the step 2i?„/v^7\^, that is decomposing each of the coordinate segments of length 2Rn 
into \/Nn parts). We will choose and fix the number iV„ later. 
For each of these parts, choose 

i_ n(/iTj(M,) 



Pj 



max 



Let ipj = 1a/j- Then, the probability that all the empirical measures /^[r„,r„+i](^j) are p^-close 
to their "theoretical" values Il{fiT„){Mj) is at least 

1 — 2ce y^ exp ' '"■' 



P,^AT„ 
16CwVaru{f,. 



As the variance Varn(^j, )(V'i) does not exceed Il{iJ,'jp^){ipj), 
probability by 



we have a lower bound for the 



N„ 



1-2ceJ2 



exp 



j=i 



iQCw ■n(/iTj(^, 



Pj 



> 1 — 2NnCE exp 



AT.. 



IQCwN^ 
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So, taking N„. = v^ ~ (AT„)3/i°, we see that the series 

n ^ " ^ n 

converges, so almost surely for all n sufficiently big, all the closeness conditions on the occu- 
pation measures are satisfied: the measures fl[T„,T„+i]{Mj) are a.s. pj-close to n(yUr„)(Mj). 

Now, let us estimate the CT„-centered distance Tp^"{fi[T„,T„+i],^{fJ'T„)), provided that these 
conditions are fulfilled. Indeed, first transport inside each Mj the part min(/i[T„,T„+i], n(/ir„)): 

we pay at most P(3L„)£:„ = O i (AT„)~5d j. Next, bring the exterior part of Il{fj.Tn) to the 

ball Uji^{cT„): due to the exponential decrease estimates, we pay at most 

as Rn = 31ogT„. Finally, let us re-distribute the parts left: we pay at most 



±p,RnP{Rn) = i?.i'(i?.)Emax(i,,H(/^|li^) 

.7=1 j=l \ n n / 

< RnP{Rn)J2\W^ + 

,=1 V^^ri 



""" ' ^ n(/iTj( M 

Nn 
< 2RnP{Rn)^ = O {{ATn)-'/') . 

Adding these three estimates, we obtain the desired Tp^" {p.[T„,T„+i],^{fJ'T„)) = O ((AT„)~'') 
with/3 = min(8Cv^, (5d)-i). ' D 

1 

Putting Lemmas [TT] and [T^ together, and recalling that AT„ ~ Tn , we conclude that almost 

surely, for all n sufficiently big, 

r;^"(/^[T„,T„,.],n(;.^j) < r„-'"^°(^^'^'^). 

Proposition [7] is proven. 

3.1.2. Euler method error control. 

Proof of Proposition 0. We prove the proposition by induction on j . The case j = i is obvi- 
ous: the only term in the right-hand side is Cih, being an estimate for the distance to the 
smoothened convolution: 

Tp^' if^f^ , ^^T,) < / / H •-P((max(|x-crJ,|x + t;-CTj)) , d/i(x) 

jR-i Juuio) vol{Uh{0)) 

< [ Pi\x-CT,\ + h)-hdfiix)<Pih)\\fi{- + CTMp-h = Cih, 
provided that h < 1 (because the norm ||/i^. ||p is bounded due to the exponential tails of fi'^). 
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Let us now check the step of induction. Namely, assume that the conclusion holds for some 
j > i, and check it for j + 1. To do this, first shift the center of the translation distance from 
ct+1 to ct '■ from Proposition E] 

Tp"^-"' (-, ■)<(! + const |cT,^i - CT,\)Tp"^ (-, ■), 
provided that |cr+i — ct\ < 1- On the other hand, we have by Lemma [6] 

|cr,+i - ct^ I < Lip^^ ^(c) ■ Tp^' (/iT,+i, /xt,) < const — ^, 

J-j+i 



ATA _CT. , . /T, 



%+i f. A^ h ^ ...of .f^ -r^v. .^ ^ fl±i -r'^: 



A 



SO finally 



(49) r/^-(-,-)< (^1 + const -^jTp^^ (-,■)< (^^J r/\-,-;. 

Now, the map 11 is Lipschitz on Ka^c by Proposition HI so for any two measures z/i, z/2 one 
has 

Tp{<!>^^\u,), $f ^(z.^)) < (i + ^(Lip^^ ^(n) + 1)) Tp{u„ u,) < {^^ ' rp(z/i, u,). 

Substituting for z/i and z/2 respectively the translated by ct^ images of measures $^(yUj.. ) 
and ht respectively, we see that 

(50) r;^^($r(4?),$r(/^^.)) < (^)^Vr^($K4?),/^T,). 

Now, using that by Proposition [71 

rr^($r(/^T,),/.T,,j<(^)(AT,)-^ 

with /3 = min(8Cvi/, (Sc?)"^), we see that 

(51) rr-H$r(/^i?^/^T,.j < (^)^Vr^($r(4?),/^T,,j < 

< {^y (7^^($r(/^s^^r(/^T,))+7^^($r(/.^),/.^,j) < 






^1+^2 /rp \ Al 



(rp \ Al+A-l / rp \ ^1 A T^ 



ii'AT-.^-/3 
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Finally, we fix the choice oi A := Ai + A2, and, using the induction assumption, the right- 
hand side of (IHTj) is not greater than 






£|^'^-'>"(1^) M^ 



The induction step is proved. 



n 



3.1.3. Decrease of energy. This section is devoted to the proof of Proposition [91 To estimate 
the decrease of energy, we will need the following 

Lemma 13. For any jj, G Ka^c, we have ffj,{fi) — V9^(n(yu)) > g{J'{fi\poo)), where 

r ^7^^, 0<E<eo<l 
9{E) = I fyieo), eo<E<ei 

( E + {g{ei)-ei), E > e, 

is an increasing continuous function, and the constants Ct^Eq^Ei depend only on a and C . 

We postpone its proof, but we use it as a motivation for the next result, which immediately 
implies Proposition [HI 

Lemma 14. There exists uq such that for any fi G Ka^c o,nd for any j > i > no: 

-F($K/^)|Poo)<l/(T,), 

where y is the solution to 

(52) y^-Y-f, 

with the initial condition y{Ti) = max(J-'(/ij;), 1). 

Proposition M is its immediate corollary, as the solution of f lS^ decreases exponentially 



for big energies y and has the form y{t) = exp{— '''^\/^{k + l)log(t/To)} for y < Eq (what 

happens for t large enough). 

We will need the following corollary to Lemma (H) 

Corollary 1. For any fixed a, C, there exists C" such that for all fi G Ka^c> for all < X < 1 , 

^((1 - x)fi + xu{fi)\p^) < ^(^Ipoo) - A ivM - </'/.(n(/i))) + c"x^ 

Proof. For fi G Ka^c, the integral that is the coefficient before A^ is uniformly bounded. D 
Let us now prove the previous lemmas. 
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Proof of Lemma\l^ Recall that, due to Corollary [H we have once /i G -ft'a.c, 

^((1 - A)/i + An(/i) \p^) < ^(^Ipoo) - A {^,{^1) - <^^(n(p))) + c"\\ 

Now note that, if hq is chosen sufficiently big, we have for any j: 

(53) C"^ < ^MM. 

Tj+i 3 

Indeed, the left-hand side of fj53|) decreases as 4, while its right-hand side decreases as 

exp{— ''+-^/^(fc + 1) logTj} ^ i. Now, for every fij := $i(/i), we have fij G -K'q.c due to 
Lemma [9] and hence due to Lemma [131 

<^/i,(/ii) - <^/i,(n(/ij)) > ^(-F(/ij|poo))- 

Hence, proving the statement of the lemma by induction on j, we have to deduce from 
•^{i^j\Poo) < yiTj) the analogous statement for /ij+i, given that 

AT- / AT\^ 1 AT- 

Hn+i\Poo) < y{T,) - g{y{T,))—L + C" -^ < y{T,) - -9{y{T,))-^. 



^ 4 Ar, 
uq to be chosen sufficiently big, we have 



Let 9j = logTj. Then, A9j := Oj^i — 9j < | j^-^ for all j large enough. So, once again asking 



Hn+iIPoo) < y{T,) - - ■ -,9{y{T,))M, = m) - '-^^^^0,, 



where y{9) = y(e^). We conclude by noticing that g{y) is an increasing function of y. So, as 
y{9) is solution to the equation y{9)' = — ^^^ , we have 

m) - ^A», < «»,«). 

hence J-'(/ij+i|poo) < vi^j+i) = viTj+i), thus proving the induction step. D 

Proof of Lemma UM Note ffist that, for p G -K"a,C; the integral Jj\fi — Il{fi)){dx)W{x — y){fi — 
n(/i))(d|/) is bounded by a uniform constant. Thus, due to LemmalU ip^{p) — ip^{Il{^)) admits 
a lower bound 

(54) ifM - Vf^i^ip)) > J^{p\poo) - Ca 

with the constant Ca being uniform over all p G Ka^c- 

Now, let us give another way to estimate the difference ^^{jj) — V9^(n(p)). Indeed, n(p) is 
the global minimiser of J-", hence for any measure p, we have 

(55) <^^(p) - <^^(n(p)) > v9^(p) - (^^(p). 

Recall that the free energy functional J-" is displacement convex. Denote by ^^ = (1 — s)^o + 
■5^) < s < 1, the quadratic Wasserstein optimal transport between p = {law of ^o} and 
Poo(' + p) = {law of i^i} and let Vg = {law of ^g}. Then, 

^(z^sIPoo) > (1 - s)J'(pIpoo). 
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Thus, we have due to Lemma [HI 

> sJ^{^\poo) + 2 ^^"^ fi){dx)W{x - y){vs - p){.dy). 

Let us now estimate the second term in the right-hand side of this inequahty. Indeed, let 
ijlo^Vi) be an independent copy of (^O)'Ci)- Then 



W{x - y){us - /i)(dx)(z/, - ^i){dy) = 

= E [Wi^o - Vo) - W{^, - r^o) - W{^o - Vs) + W{^, - r/,)] . 

For any fixed L, we can divide this expectation into two parts: the one corresponding to 
max (|^j|, |?7j|) > L and the one with |^j|, \rji\ < Lioi i = 0, 1. We also remind that z/j e -f^o,C2 

ij6{0,l} 

for i = 0, 1 and that P controls W as well as its first and second derivatives. So, there exists 
a positive constant C such that 

|E [W{^o - Vo) - W{^s - Vo) - W{^o - Vs) + W{^s - Vs)] I 



/oo 
H^(2/)dF^ax«„,^,,,„,,,)(/) 



max (|g,|,|r;j|)<L 
»,J6{0,1} 



max P( la; 1)1^0 — CsWvo — Vs\^ 

|x|<4L ^' '^'^" ^''"" '"' \ max {|?,|,|r,,|)<L 



< E 



oo 



+ A I P{2l)d{l-C2e 



So, using the already mentioned comparison W2{fi,Poo) < (f~-^(/^|Poo), we have 



,2 D//1 rMl/2/,. „ \ r<TDlnT\„-OLL 



VM - V^/.(^.) > 5-^(/i|Poo) - s'P{AL)Wi{^i, p^) - CP{2L)e- 



We decide from now on to fix s = jpu^s^ with the choice of L to be fixed later. Then, 
s - ^s'^P{AL) = s/2 and 

(56) ipM - M^s) > ^^fif'lPoo) - CP{2L)e--^. 
For J^(/i|poo) sufficiently small, fixing L = ^\ log J-'(/i|poo)|, we have 

(57) 16^'^^^'''°°) - ^^(2^)^""' 
and hence the right-hand side of (!56|) is estimated from below by 

^ -^(pIPoo) > |. ^-^/^ m. -^(pIPoo)- 



16P(4L) ^^'^°°^- |logJ^(p|p 



oo / 
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So taking g{E) := .^^'^^.^. E for such values oi E = J'{fi\pao), we have for such £"s the conclusion 
of lemma satisfied. Next, fixing eo to be such that ( 1571) is satisfied for J^(yu|poo) < £^0, and for 
any J^(/i|poo) > £o, choosing the same L as for J^(/i|poo) = ^o, we have 



what allows to deduce 

(58) c,,(p) - ^,m,)) > I ^(^^(,^) if ^(^1^^) > ,^ 

Finally, taking the maximum between the right-hand side of (l5lll and (l58l) . we obtain the 
desired conclusion. D 

3.2. Proof of Theorem |4l. As it has been already shown in ( !22l) . we have 

Jt„ '-'Wt JTn '^Wt J-n 

Thus, almost surely one has oscti^[T,^^Tnj^i\Ct — )• as n ^^ oo. So, to prove Theorem HJ it suffices 
to show that the sequence ct„ converges almost surely. 
Now, let us estimate the distance ct„+i — ct^- Indeed, 

/^T„+i = /iT„ + 7^ — (/i[T„,r„+i] - /^r„)- 

Translating ct„ to the origin, using the decrease estimates of §3.11 and recalling that c(-) : 
K^ ^ — )■ M*^ is a 7p-Lipschitz function, we see that 

|ct„+i - ctJ < Lip^o (c) ■ — ^ ■ 7^^"(/i[T„,r„+i],/iTj. 
As in §3.1.2[ the distance in the right-hand side can be estimated as a sum of two distances: 

(59) r;"" (/i[T„,T„+,], /iTj < r;"" (/U[t„,t„+i], n(/iT^j) + r;-- (n(/iTj, ptJ. 

We already have an estimate for the first term in this sum: 

r;-(/^[T„,T„,,,n(/.))<T„-'^^'^^^^-'^). 

On the other hand, the limit density poo is a fixed point of the map 11. And the map 11 being 
Lipschitz on K^(j, the second summand in fl59l) can be estimated as 

r;""(n(pTj,/iTj < (Lipxo ^(n) + 1) ■ rp(p^„,poo). 

The latter distance is already estimated in the proof of Theorem [3l almost surely for n 
sufficiently big, we have 

rp(p^„,Poo) < exp{-a ^Vlog^n}. 
Finally, adding the estimates for the first and the second terms in (1591) . we obtain that for all 
n sufficiently big, 

r;""(/^[T„,T„,,],/XTj < Tn'^'-^^'''''^^ + (Lip^o jn) + l)exp{-a ^Vi^^} 
and hence 
|CT„,, - ctJ < Lip^o Jc) ■ ^ /'t:-'^^^-'^) + (Lip^o^(n) + l)exp{-a '^^/\^^n]] • 
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We choose T„ = r?!"^ and so ^^ x n ^ Hence 



1 exp{-a'=+^y|logn} 
2^ |cT„+i - ctJ < const }_^ ^i+^i„(4c^,i/(io^)) + const ^ 



n 



Both series in the right-hand side converge, and thus the series ^„ |ct„+i — ct„| converges 
almost surely. This concludes the proof. 



Appendix 1: Singularity at t = 

Let us now prove that a solution to the equation ([1]) with any initial condition at t = 
(where the equation has a singularity) exists and is unique. 

Proposition 11. For any xq and almost every trajectory Bt of the Brownian motion, a (con- 
tinuous at t = 0) solution Xf to the equation (^ with the initial condition Xq = xq exists on 
all the interval [0, +oo) and is unique. 

Proof. As Proposition [5] provides us global existence and uniqueness of solutions, starting from 
any arbitrary positive time r > 0, it suffices to check the existence and uniqueness on some 
interval [0, 6). For the sake of simplicity of notation, suppose that xq = 0. 

Let 6i > Ohe such that for all < t < 6i, \Bt\ < \ and 5i sup|2.|<2 |Viy(a;)| < |. We work 
on the trajectories, which are staying inside Wi(0), the unit ball centered in xq = 0. So, we 
consider X, : [0, 6i) — )■ Wi(0), t (-)■ Xt. Denote by fif the empirical measure of the process X. 
Then, the application x : X i— )■ A is such that 



X = R 



VW*fif{X,,)ds 



is well-defined on this space, and Xt also remains stuck in Ui{0). Indeed, for any time t < 5i, 
such that the solution A is defined on [0,t] and stays in Wi(0), we have 



VW*ijf{Xs)ds 



- f VW{X,-Xu)duds 
'0 -5 Jo 

/■* 1 /■' 1 

< / -/ snp\VW{x)\duds<6isnp\VW{x)\<-^ 

Jo S Jo |x|<2 |x|<2 3 



Thus, if there existed a time to ^ "^i such that {Xt^l > 7/8 for the first time, then we would 
see that \Xtg\ < 1/2 + 1/3, which would contradict the bound lAi^l > 7/8. So, A stays in 
Ui{0) for any <t<6i. 

Let us now show that for 6 < 6i sufficiently small, the map x is a contraction on the space 
of continuous maps A, from [0, 6] to Wi(0) with Aq = 0. Indeed, consider now two trajectories 
A*^^) and X^'^\ realizing a coupling with the same Brownian motion, and their respective 
images (by x) ^^^^ and X^'^\ Then, denoting by Lip(iy) the Lipschitz constant of VW on 



ERGODICITY OF SELF-ATTRACTING MOTION 33 

2t 



the 2-radius ball, Lip(iy) := sup{||VW(x)|| : |x| < 2}, we have 



|x«-xf)| 



< 



Jo Jo 

[ - r VW^(Xii) - X^i)) - VW{Xi^^ - X^2)) duds 
Jo -s Jo 

[ - r \VW{Xi^^ - X^i)) - VW{Xi^^ - Xi^^)\du ds 

Jo -5 JO 

f - /\ip(W^)(|Xi^) - XP\ + |X(i) - X(^)|) duds 
Jo ^ Jo 



< 

< tLip(W^)(||X« -X(2)||c([o,5]) + ||X« -X(2)||^(fo^,])), 
where ||X||c([o,5]) is the norm of X on the space C([0,(5]). As t < 6, we conclude that 
||X(i) - X(2)||^([o_,]) < 6 Lip(iy)(||XW - X(2)||^(p,,]) + ||X(i) - X('^\\ci[o,5]))- As soon as 
6 Lip(W^) < 1, we have 

We choose S such that (5 Lip(iy) < 1/3 and then x is a contraction, as stated, with Lip(x) < 
1/2. So, we have obtained existence and uniqueness of the solution on [0, 6]. D 

Appendix 2: Non-symmetric counter-example 

We end this paper with an example showing that for a non-symmetric interaction potential 
W, the conclusion of Theorem [1] does not hold. 

Consider a non-symmetric quadratic interaction potential W{x) = |(x — 1)^. Then, the 
averages of the process {Xt)t defined by ([1]), j J^ Xsds = q — 1 tend to -foo. 

To motivate this behaviour, heuristically, we first note that, for any finite- variance measure 
u, the convolution W * u equals 

1 1 x^ 
W * iy{x) = -(x - 1)2 - (x - l)E(z/) + -E(z/2) = (E(z/) + l)x + const 

and hence n(z/) is the Gaussian law A/'(l + E(z/), 1). Thus, if we consider a trajectory of the 
approximating flow Ut = j(n(z/t) — z/j), we have for its mean value 

dt t^ ' V 

and so Eut ~ logt. 

For a formal proof, note that (as the interaction potential is a polynomial of degree 2) the 
evolution of the couple (Xj, q), where q = c{nt) = E/i^ + 1 is Markovian: 

dXt = dBt - {Xt - ct)dt, 

Q = i(X,-Q + l). 

Changing Xf to Yf = Xf — Ct, we obtain 

dYt = dBt-{Yt + \{Yt + l))dt, 

Ct = \{Yt + l). 

The equation on Y does not contain q. So, explicit solution of this system and rigorous 
justification of the desired properties become an exercise. 
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